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Abstract. We study the angular bispectrum of local type arising from the (possibly cor- 
related) combination of a primordial adiabatic mode with an isocurvature one. Generically, 
this bispectrum can be decomposed into six elementary bispectra. We estimate how pre- 
cisely CMB data, including polarization, can enable us to measure or constrain the six cor- 
responding amplitudes, considering separately the four types of isocurvature modes (CDM, 
baryon, neutrino density, neutrino velocity). Finally, we discuss how the model-independent 
constraints on the bispectrum can be combined to get constraints on the parameters of 
multiple- field inflation models. 
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1 Introduction 

Inflation is currently the best candidate to explain the generation of primordial perturbations, 
but many of its realizations remain compatible with the present data. One can hope that 
future data will enable us to find additional information in the primordial perturbations that 
could help to discriminate between the various mechanisms that can have taken place in the 
very early Universe. 

In this respect, it is important to test the adiabatic nature of the primordial pertur- 
bations. Since single-field inflation predicts only adiabatic perturbations, the detection of a 
fraction of an isocurvature mode in the cosmological data would rule out the simplest mod- 
els of inflation. By contrast, multiple-field inflation could easily account for the presence of 
isocurvature modes [1], which can even be correlated with the adiabatic component [2, 3]. 

As shown in [4] , the most general primordial perturbation is a priori a linear combination 
of the usual adiabatic mode with four types of isocurvature modes, respectively the Cold 
Dark Matter (CDM), baryon, neutrino density and neutrino velocity isocurvature modes. 
The existence, and amplitude, of isocurvature modes depends on the details of the thermal 
history of the Universe. Various scenarios that can lead to observable isocurvature modes 
have been discussed in the literature (double inflation [2, 5, 6], axions [7], curvatons [8-12]). 

In parallel to the possible presence of isocurvature modes, another property that could 
distinguish multiple-field models from single-field models is a detectable primordial non- 
Gaussianity of the local type. So far, the WMAP measurements of the CMB anisotropics [13] 
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have set the present limit f^ al = 32±21 (68 % CL) [and -10 < f^ 1 < 74 (95 % CL)] on the 
parameter /nl^ that characterizes the amplitude of the simplest type of non-Gaussianity, 
namely the local shape. Similarly to isocurvature modes, a detection of local primordial 
non-Gaussianity would rule out all inflation models based on a single scalar field, since they 
generate only unobservably small local non-Gaussianities [14]. Scenarios with additional 
scalar fields, such as another inflaton (see e.g. [15, 16]), a curvaton [17] or a modulaton [18- 
20], which can produce detectable local non-Gaussianity, would then move to the front stage. 

Isocurvature modes are usually investigated by constraining the power spectrum of 
primordial perturbations with CMB or large-scale structure data (see e.g. [21-28]). How- 
ever, isocurvature modes could also contribute to non-Gaussianities as discussed in several 
works [29-37]. Moreover, there exist models [35] where isocurvature modes, while remaining 
a small fraction at the linear level, would dominate the non-Gaussianity. As shown in [37], 
these CDM isocurvature modes would be potentially detectable via their non-Gaussianity in 
the CMB data such as collected by Planck. Non-Gaussianity can thus be considered as a 
complementary probe of isocurvature modes. 

In the present work, we refine and extend our previous analysis [37] by considering all 
types of isocurvature modes, not only the CDM isocurvature mode. We analyse the bispec- 
trum generated by the adiabatic mode together with one of the four isocurvature modes. The 
total angular bispectrum can be decomposed into six distinct components: the usual purely 
adiabatic bispectrum, a purely isocurvature bispectrum, and four other bispectra that arise 
from the possible correlations between the adiabatic and isocurvature mode. Because these 
six bispectra have different shapes in /-space, their amplitude can in principle be measured 
in the CMB data and we have computed, for each type of isocurvature mode, the associated 
6x6 Fisher matrix to estimate what precision on these six parameters could be reached with 
the Planck data. We also show that the inclusion of polarization measurements improves the 
predicted precision of some isocurvature non-Gaussianity parameters significantly. 

The elementary bispectra discussed above depend only on the adiabatic and isocurvature 
transfer functions and are thus independent of the details of the generation mechanism. Now, 
by assuming a specific class of inflationary models, one obtains particular relations between 
the six bispectra, which can be used as consistency relations for the model or to constrain 
the model parameters. We illustrate this in the context of curvaton-type models, generating 
adiabatic and CDM isocurvature perturbations. 

The outline of the paper is the following. In the next section, we present the various 
isocurvature perturbations and discuss their impact on the CMB angular power spectrum. 
The following section is devoted to the angular bispectrum and its decomposition into six 
elementary bispectra. We then discuss the observational prospects to detect these elementary 
local bispectra in the future data, distinguishing the various isocurvature modes. Finally, we 
consider models where primordial perturbations are generated by an inflaton and a curvaton, 
and show that the amplitudes of all six bispectra depend on only two coefficients, which can 
be constrained from the data. The last section contains our conclusions. 

2 Isocurvature perturbations 

In this section, we recall the definition of isocurvature modes in the context of linear cosmo- 
logical perturbations. At the time of last scattering, the main components in the Universe are 
the CDM (c), the baryons (b), the photons (7) and the neutrinos (v). All these components 
are characterized by their individual energy density perturbation 5pt and their velocities Vj 
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(as well as higher momenta of their phase space distribution functions, which we do not dis- 
cuss here; see [38] for details). The "primordial" perturbations for each Fourier mode k are 
usually defined on super-Hubble scales, i.e. when k <C aH, deep in the radiation dominated 
era. 

The most common type of perturbation is the adiabatic mode, characterized by 

5n c _ 5n b _ 5n v _ <5n 7 

— — — , v^-J-j 

n c rib n v n ~j 

which means that the number of photons (or neutrinos, or CDM particles) per baryon is 
not fluctuating. In terms of the energy density contrasts (5 = 5p/p), the above condition is 
expressed as 

3 3 

5 C = 5b = -5 U = -<5 7 , (2.2) 

where the 3/4 factor, for photons and neutrinos, comes from the relation p oc ra 4 / 3 for 
relativistic species. 

Assuming adiabatic initial conditions is natural if all particles have been created by the 
decay of a single degree of freedom, such as a single inflaton, and, so far, the CMB data are 
fully compatible with purely adiabatic perturbations. However, other types of perturbations 
can be included in a more general framework. In addition to the adiabatic mode, one can 
consider four distinct types of so-called isocurvature modes [4]: the CDM isocurvature mode 
S c , the baryon isocurvature mode Sb, the neutrino density isocurvature mode S u d and the 
neutrino velocity isocurvature mode S vv . At zeroth order in kr, where r is the conformal 
time, the first three are characterized, respectively, by 





3 A 


3 A 


5 b = 




= 4^7 


5 C = 






4*" 


= 4*7 
3 A 


5 b = 


= 6 C -- 


= 4*7 



S c = S c + -j&y, 5b = ~ A 5 V = -5 7 (CDM isocurvature) (2-3) 
5b = Sb + -<5 7 , 5 C = -5 U = ~6ry (baryon isocurvature) (2-4) 
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3 3 

-5 U = S u d + -5 7 , 5b = 5 C = -5y (neutrino density isocurvature), (2.5) 

while the corresponding velocities tend to zero. As for the neutrino velocity isocurvature 
mode, it is characterized by non vanishing "initial velocities", 

7 /4\ 4/3 

V y = S vv , V 7 6 = —-N u I — S uv (neutrino velocity isocurvature), (2.6) 

8 V 11 / 

where V 7 b is the common velocity of the photon-baryon plasma (the photons and baryons are 
initially tightly coupled via the Thomson scattering off free electons) and N v is the number 
of species of massless neutrinos. The above relation between the two velocities ensures that 
they exactly cancel in the momentum density, while the energy densities satisfy the adiabatic 
condition (2.2). 

In the following, each mode will be characterized by its amplitude: the curvature per- 
turbation on constant energy hyper surf aces, £, for the adiabatic mode, and S c , Sb, S u d and 
Suv f° r the four isocurvature modes. These variables will be denoted collectively as X 1 . In 
the context of inflation, a necessary condition for these isocurvature modes to be created is 
that several light degrees of freedom exist during inflation. Moreover, since the adiabatic 
and isocurvature modes can be related in various ways to these degrees of freedom during 
inflation, one can envisage the existence of correlations between these modes. 
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These various modes lead to different predictions for the CMB temperature and polar- 
ization. Let us consider for instance the temperature anisotropics, which can be decomposed 
into spherical harmonics: 

-jT = S ^ai m Yi m . (2.7) 

lm 

The multipole coefficients a\ m can be related linearly to any of the primordial modes. The 
precise correspondance can be computed numerically and written in the form 

a\ m = 47TH)' / (0T*'( k ) dlik) ^(k) , (2.8) 

where g\ (k) is the transfer function associated with the corresponding primordial perturba- 
tion {gj{k) depends also on the various cosmological parameters). 

For each type of perturbation, the angular power spectrum is thus given by 

9 poo 

C{ = (aLO = dkk 2 [gfikfMk), (2.9) 

K Jo 

where we have introduced the primordial power spectrum Pi(k) defined by 

(X^X 1 ^)) ee (2vr) 3 5(k 1 + k 2 )P / (fc 1 ) . (2.10) 

For our purposes, the crucial point is that the transfer functions associated with isocurva- 
ture perturbations are very different from the adiabatic transfer function. Moreover, each 
isocurvature mode leads to a specific signature that enables one to distinguish it from the 
other isocurvature modes. The only exception are the CDM and baryon isocurvature modes 
which give exactly the same pattern, up to the rescaling: 

Sb = UbcS c , LObc = ^ , (2-11) 

where the parameters Qb and f2 c denote, as usual, the present energy density fractions, 
respectively for baryons and CDM (note however that these two modes can in principle be 
discriminated via other effects, see e.g. [25, 39-41]). 

More generally, when we also allow for possible correlations between the modes and 
include E-polarization, the angular power spectra are given by 

c lJ;aa' = i dk k 2 gj'* (k) g^' (k) Pj J (k) , (2.12) 

n Jo 

where I and J label the isocurvature mode and a and a' the polarization (i.e. either T 
(temperature) or E (polarization)). The primordial power spectrum Pjj(k) is now defined 
by 

(X / (k 1 )X J (k 2 )) ee (27r) 3 ,5(k 1 + k 2 )Pu(h) , (2.13) 

which generalizes the previous definition to include the presence of correlations between the 
modes, which corresponds to the situation where Ppj with different / and J does not vanish. 

All this is illustrated in Fig. 1 and 2, where we have plotted the angular power spectra 
for all the various modes (Fig. 1) and the isocurvature cross power spectra where one of the 
components is adiabatic (Fig. 2), assuming the same primordial power spectrum for all. As 
one can see from these figures, the CDM (and baryon) isocurvature mode decreases much 
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Figure 1. Angular power spectra (multiplied by 1(1 + 1)) for the temperature (left) and polarization 
(right) obtained from purely adiabatic or purely isocurvature initial conditions. The amplitude and 
spectral index of the primordial power spectrum, as well as the cosmological parameters, on which 
the transfer functions depend, correspond to the WMAP7-only best-fit parameters (WMAP7-only 
best-fit parameters are used for all the figures and explicit computations in this paper). 



cross power spectra 




EE cross power spectra 




Figure 2. The isocurvature cross power spectra (multiplied by 1(1 + 1) and in absolute value) with 
one component fixed to be the adiabatic one. The figure on the left shows temperature (TT), the one 
on the right polarization (EE). 



faster with I than the other modes. In fact it turns out that if one multiplies the CDM 
isocurvature power spectrum by I 2 (I + l) 2 instead of 1(1 + 1), it falls off roughly in the same 
way as the other modes at large I, as illustrated in Fig. 3. This figure also nicely shows the 
relative phases of the acoustic peaks for the different modes. 

When the "primordial" perturbation is a superposition of several modes, the multipole 
coefficients depend on a linear combination of the "primordial" modes, 

a lm = ^(-i) l J ^ ^£*'(k)rf(fc)J Yf m (k). (2.14) 

(Here we have once again omitted the polarization indices, as we will do in most of the 
equations of the paper, in order to improve readability.) As a result, the total angular power 
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Figure 3. Same angular power spectra for the temperature and the polarization as in Fig. 1, except 
that the horizontal axis is now linear, and that the CDM isocurvature spectrum is multiplied by 
10~ 4 Z 2 (7 + l) 2 instead of 1(1 + 1). 



spectrum is now given by 

2 f°° 

d = {a lm a* lm ) = V - / dkk 2 g{(k)gf(k)P I j(k). (2.15) 

We infer from CMB observations that the "primordial" perturbation is mainly of the adiabatic 
type. However, this does not preclude the presence, in addition to the adiabatic mode, of 
an isocurvature component, with a smaller amplitude. Precise measurement of the CMB 
fluctuations could lead to a detection of such an extra component, or at least put constraints 
on its amplitude. For example, constraints on the CDM isocurvature to adiabatic ratio, 

a = -f, (2.16) 

based on the WMAP7+BAO+SN data, have been published for the uncorrelated and fully 
correlated cases (the impact of isocurvature perturbations on the observable power spec- 
trum indeed depends on the correlation between adiabatic and isocurvature perturbations, 
as illustrated in [3]). In terms of the parameter a = a /(I + a), the limits given in [13] are 1 

a < 0.064 (95%CL), a x < 0.0037 (95%CL) , (2.17) 

respectively for the uncorrelated case and for the fully correlated case. 



3 Generalized angular bispectra 

In this section, we turn to non-Gaussianities, including both adiabatic and isocurvature 
modes. 

x Our notation differs from that of [13]: our a corresponds to their a and our fully correlated limit corre- 
sponds to their fully anti-correlated limit, because their definition of the correlation has the opposite sign (see 
also [42] for a more detailed discussion) . 
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3.1 Reduced and angular- averaged bispectra 

The angular bispectrum corresponds to the three-point function of the multipole coefficients: 

r>m,im2rii3 i • \ /o i \ 

hhh = \ a hm 1 a l2m2 a l3m3 / • K^- 1 ) 

Substituting the expression (2.14) into the angular bispectrum, one can write it in the form 

Kl7r 3 = Sr^ya , (3-2) 
where the first, purely geometrical, factor is the Gaunt integral 

<cr m3 = / d2fi *Wi y ^ ^ 3 (n) . (3.3) 

The second factor, usually called the reduced bispectrum, is given by 

l-OO 

J2 



r drj h (kir)j l2 (k 2 r)j l:i (k 3 r), (3.4) 

which depends on the bispectra of the primordial X : 

(X I (k 1 )X J (k 2 )X K (k 3 )) = (27r) 3 d('S i k i )B IJK (k 1 , k 2 , k 3 ) . (3.5) 

The reduced bispectrum (3.4) is the sum of several contributions, corresponding to different 
values of the indices /, J and K that vary over the range of modes included in the primordial 
perturbations. This expression thus generalizes the purely adiabatic expression given in [43] . 
It is also useful to define the angle- averaged bispectrum 



rj = X- fh h h\ umi m 2m3 _ / (2Zi + l)(2i 2 + l)(2Z 3 + l) (hhk\ h 

*hhh - 2^ { mi m 2 m 3 ^Wa'a ~ V 4tt I J bhhh ' 

(3.6) 

where the second relation is obtained by substituting (3.2) and by using the identity 



V I h h h \ cmi m 2 m 3 = / (2jl + 1) (2Z 2 + 1) (2l 3 + 1) / h h h \ 
^— ' \ nriA mn m,n 



mi m 2 m 3 J ' l(2 ' 3 V 4tt \ 

mi,m2,m3 x ' v 

3.2 Non-Gaussianities of local type 

To proceed further, one must make some assumption about the functional dependence of the 
bispectra B IJK (k\, k 2 , k 3 ) in Fourier space. This corresponds to the so-called "shape" of the 
bispectrum [44], which has been discussed at length in the literature in the purely adiabatic 
case where the B IJK reduce to the single bispectrum B^. In the present work, we consider 
the simplest form of non-Gaussianity, namely the local shape. In the purely adiabatic case, 
it is defined by 

C(x) = Cl(x) + ^/ NL (Cl(x) 2 - (Cl> 2 ) , (3.8) 

in physical space, where the factor 3/5 appears because /nl was originally defined with 
respect to the gravitational potential $ = (3/5)C, instead of £. The subscript L here denotes 
the linear part of the perturbation, which is assumed to be Gaussian. 
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In Fourier space, this leads to a bispectrum that depends quadratically on the power 
spectrum: 

£ ?CC = /nl [P c (k 2 )P c (k 3 ) + P c (k 3 )P c (ki) + P c (*i)P c (*2)] , /nl = ^/nl • (3.9) 

In the present context where we assume the presence of an isocurvature mode in addition 
to the dominant adiabatic mode, the simplest extension of (3.9) is to assume that all the 
generalized bispectra B IJK (ki,k 2 ,k 3 ) can be written as the sum of terms quadratic in the 
adiabatic power spectrum (note that this implicitly assumes that the power spectrum of the 
isocurvature mode and the isocurvature cross power spectrum, if nonvanishing, have the same 
spectral dependence as the adiabatic one). However, in contrast with (3.9) where all terms 
share the same coefficient, as a consequence of the invariance of the bispectrum under the 
exchange of momenta, this is no longer the case for the generalized bispectra when the indices 
/, J and K are not identical. What the definition (3.5) implies is simply that the bispectra are 
left unchanged under the simultaneous change of two indices and the corresponding momenta 
(e.g. I and J, k\ and k 2 ). This leads to the decomposition 

B IJK (h, k 2 , k 3 ) = f^ K P ( (k 2 )P ( (k 3 ) + fJf'P^Pdh) + f«i IJ P c (h)Pdk2) , (3.10) 

where the coefficients /J' L must satisfy the condition 

f^ K = f^L J - (3-11) 

To keep track of this symmetry, we separate the first index from the last two indices with a 
comma. 

3.3 Link with multiple-field inflation 

It is instructive to show that our definition of generalized local non-Gaussianity is the natural 
outcome of a generic model of multiple-field inflation. Indeed, allowing for several light 
degrees of freedom during inflation, one can relate, in a very generic way, the "primordial" 
perturbations X 1 (defined during the standard radiation era) to the fluctuations of light 
primordial fields 4> a , generated at Hubble crossing during inflation, so that one can write, up 
to second order, 

X 1 = Nl 6^ + l -Nl b 5<f> a S<f, b + ... (3.12) 
where the d(p a can usually be treated as independent quasi- Gaussian fluctuations, i.e. 

<<^ a (k) «5^(k')) = (2vr) 3 5 ab P H (k) 5(k + k') , P S4> {k) = 2ir 2 k- 3 2 , (3.13) 

where a star denotes Hubble crossing time. The relation (3.12) is very general, and all the 
details of the inflationary model are embodied by the coefficients N[ and iV^. 

Substituting (3.12) into (3.5) and using Wick's theorem, one finds that the bispectra 
Bijk can be expressed in the form 

B IJK {k u k 2 , fc 3 ) = X'^Ps^Ps^h) + X J ' KI P s<i> {h)P S(j> (k 3 ) + \ KJJ P&^Ps^) , 

(3.14) 

with the coefficients 

x i,jk _ s ac 5 bd Nl b N c J Nf (3.15) 
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(the summation over scalar field indices a, b, c and d is implicit), which are symmetric under 
the interchange of the last two indices, by construction. Since the adiabatic power spectrum 
is given by 

P c = (5 ab N^N c b )P s ^ = AP s</> , (3.16) 

one obtains finally (3.10) with 

ftt K = 4 J L K /A 2 , (3.17) 

where it is implicitly assumed that the coefficients are weakly time dependent so that 
the scale dependence of A 2 can be neglected. One can notice that the first index is related 
to the second-order terms in the decomposition (3.12), while the last two indices come from 
the first-order terms. 

Except in the last section devoted to a specific class of early Universe models, all our 
considerations will simply follow from our assumption (3.10) and thus apply to any model 
leading to this local form, whether based on inflation or not. 

3.4 Decomposition of the angular bispectrum 

After substitution of (3.10) into (3.4), the reduced bispectrum can finally be written as 

I,J,K 

where each contribution is of the form 2 



with 



poo 

#££ = 3/ rHra\ h {r)Pl(r)pF 3) {r), (3.19) 

J 



«f(r) = ~ / *?dkMkr)gf{k), (3.20) 



7T 

Pi(r) = \\ k?dkji{kr)gfik)Pc(k). (3.21) 

While we have omitted the polarization indices, the reader should keep in mind that each 
transfer function carries, in addition to the isocurvature index, a polarization index, and 
hence the same is true for a/ and fy. As a consequence, the bispectrum has three polarization 
indices that we do not show. 

The purely adiabatic bispectrum, usually the only one considered, can be expressed as 



b' 



oo 

Ws = 3 / r 2 dra\ h (r)^(r)(3f 3) (r) 



oo 

2 



r dr 



of h {r)Pi (r)# 3 (r) + «<■ (r)tf 3 (r)/% (r) + (r)^ (r)/3< (r) . (3.22) 



The functions a\{r) and /3f(r) are plotted respectively in Fig. 4 and Fig. 5. We have consid- 
ered both the temperature and polarization transfer functions. The radial distance r can be 
expressed as the speed of light times the difference in conformal time between now and the 
time in the past we consider. In the figures we have chosen some sample values around the 
time of last scattering, which corresponds tor ~ 14100 Mpc using the WMAP7-only best-fit 
parameters. 
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Alpha adiabatic temperature 



e _ 12 Alpha adiabatic polarization 



r=14270.4 
r=14214.7 
r=14145.0 
r=14043.8 





Figure 4. The adiabatic aj (r) as a function of I for temperature (left) and polarization (right). It 
has been evaluated at four different values of r: 14043.8, 14145.0, 14214.7, and 14270.4. 



Beta adiabatic temperature 



Beta adiabatic polarization 



r=14270.4 

- - r=14214.7 
-- r=14145.0 

— r=14043.8 








r=14270.4 

■■- r=14214.7 
-- r=14145.0 
r=14043.8 




Figure 5. The adiabatic £(Z + l)/3j f (r) as a function of I for temperature (left) and polarization (right) 
(for the same values of r as in Fig. 4). 

Since we consider local non-Gaussianity, the main contribution to the bispectrum comes 
from the squeezed limit, i.e. one of the multipole numbers is much smaller than the other 
two. To simplify the analysis, let us assume that l\ <C \<i = 1% = I. One finds that, in this 
limit, the integrand in (3.22) is dominated by (twice) aj(r)pf(r)pf(r), whereas the first term 

is negligible: one can see from Fig. 4 and 5 that pj (r) is much larger (in absolute value) than 

aj (r). This is true both for temperature and polarization (one should also keep in mind that 
the /3's have been multiplied by 1(1 + 1) in the plots, which makes them look much larger at 
large I). 

The purely isocurvature bispectrum has exactly the same structure as (3.22), but with 



the functions af(r) and /3f(r) replacing oq(r) and f}\ \{r). Moreover, the shapes of (r) and 
f3f(r) depend on the type of isocurvature mode: Fig. 6 and Fig. 7 correspond to the CDM 
isocurvature mode, Fig. 8 and Fig. 9 to the neutrino density isocurvature mode and Fig. 10 
and Fig. 11 to the neutrino velocity isocurvature mode. The functions otf(r) and Pf(r) for 



We use the standard notation: (hhh) = [hfoh + 5perms]/3! 
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Alpha CDM isocurvature temperature 



Alpha CDM isocurvature polarization 



r=14270.4 
r=14214.7 
r=14145.0 
r=14043.8 





Figure 6. The CDM isocurvature af°(r) as a function of / for temperature (left) and polarization 
(right) (for the same values of r as in Fig. 4). 



Beta CDM isocurvature temperature 




Beta CDM isocurvature polarization 




Figure 7. The CDM isocurvature 1(1 + l)/3f c (r) as a function of I for temperature (left) and polar- 
ization (right) (for the same values of r as in Fig. 4). 



the baryon isocurvature mode can be deduced from the CDM isocurvature functions by a 
simple rescaling, according to (2.11): 



a? b (r) = uj bc a? c (r) , ft b (r ) = u bc /3f c (r ) 



(3.23) 



The other bispectra depend on a mixing of the adiabatic and isocurvature functions. 
For example, one finds 



til + b i£s = 6 / r2dra k ( r K ( r K) w 



r dr 



(3.24) 



+4 2 (r) [^(r)^(r) + ^(r)^(r)] +af s (r) [^(r)^ (r) + ^(r)^(r)] } . 

Since &j j j and { cannot be distinguished, we will always consider the sum of the two, 
and similarly for and . 
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le _ 11 Alpha neutrino density isocurvature temperature 



e-12 Alpha neutrino density isocurvature polarization 



r=14270.4 
r=14214.7 
r=14145.0 
r=14043.8 





Figure 8. The neutrino density isocurvature af" d (r) as a function of I for temperature (left) and 
polarization (right) (for the same values of r as in Fig. 4). 



Beta neutrino density isocurvature temperature 




Beta neutrino density isocurvature polarization 




Figure 9. The neutrino density isocurvature 1(1 + " d (r) as a function of I for temperature (left) 
and polarization (right) (for the same values of r as in Fig. 4). 



In summary, after integration over r of these various combinations of a and j3 functions, 
we obtain six independent bispectra, for each type of isocurvature mode. To illustrate the 
typical angular dependence of these bispectra, we have plotted them as functions of I3, for 
fixed values of l\ and I2, respectively in the CDM isocurvature case (Fig. 12), in the neutrino 
density isocurvature case (Fig. 13) and in the neutrino velocity case (Fig. 14). We plot only 
the pure temperature (TTT) and pure polarization (EEE) bispectra, but of course one also 
has all the polarization cross bispectra. As mentioned before, the (£, (S) curve corresponds 
to twice bj^i since we consider the sum of and bf^; , which cannot be distinguished. 

The same applies to the (S, CS) curve. As usual, the bispectra for the baryon isocurvature 
mode are deduced from the CDM bispectra by the appropriate rescalings. 

4 Observational prospects 



As we will see in the next section, one can envisage early Universe scenarios that generate sig- 
nificant isocurvature non-Gaussianity, which could dominate the purely adiabatic component, 
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Figure 10. The neutrino velocity isocurvature af "* (r) as a function of I for temperature (left) and 
polarization (right) (for the same values of r as in Fig. 4). 



Beta neutrino velocity isocurvature temperature 




Beta neutrino velocity isocurvature polarization 




Figure 11. The neutrino velocity isocurvature 1(1 + l)f3f" v (r) as a function of I for temperature (left) 
and polarization (right) (for the same values of r as in Fig. 4). 



even if the adiabatic mode is dominant in the power spectrum as required by observations. 
This is why it is important to assess how precisely one can hope to measure and to discrim- 
inate the various isocurvature bispectra in the future. 

The most general analysis would require to consider simultaneously all five possible 
modes, which corresponds to a total of n 2 (n + l)/2 = 75 coefficients, taking into account the 
symmetry (3.11). In order to simplify our analysis, we will consider separately the various 
isocurvature perturbations. In other words, we will assume that the primordial perturbation 
is the superposition of a dominant adiabatic mode and of a single isocurvature mode. In this 
case, the total bispectrum is characterized by six parameters, which we now denote /w, 



E/ 



(4.1) 



where the index i varies between 1 to 6, following the order indicated in the upper line. 
Note that, because of the factor 2 in front of f^'^ s and f S '^ s , we define b\ 2 ) , = 26f'f^ and 

' 'l'2'3 '11213 
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Figure 12. Plot of Ify'j j | in the CDM isocurvature case, as functions of £3, for h = h = 1300. 
The figure on the left shows the temperature-only TTT bispectra, the one on the right the pure 
polarization EEE bispectra. 




Figure 13. Plot of l&f'ff^l in the neutrino density isocurvature functions of Z 3 , for h = h = 

1300. The figure on the left shows the temperature-only TTT bispectra, the one on the right the pure 
polarization EEE bispectra. 



b) , , = 2b, , , whereas there is no such factor 2 for the other terms. 
4.1 The Fisher matrix 

To estimate these six parameters, given some data set, the usual procedure is to minimize 

X 2 = (^(B obs - ^ f {i) B^), (B obs - Y f {i) B®)\ . (4.2) 

For an ideal experiment (no noise and no effects due to the beam size) without polarization, 
the scalar product is defined by 

(B,B')= Bhhl f' hhh . (4.3) 
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Figure 14. Plot of Ibj^^J in the neutrino velocity isocurvature case, as functions of I3, for h = I2 = 
1300. The figure on the left shows the temperature-only TTT bispectra, the one on the right the pure 
polarization EEE bispectra. 



The bispectrum variance in that case is given by 

a hhh = ( B hhh) ~ ( B hhh) 2 ~ ^hhh C h C h C h ( 4 - 4 ) 
in the approximation of weak non-Gaussianity, where 

A llhh = 1 + S hh + S hh + + 2 S hh S hl3 ■ ( 4 - 5 ) 

The best estimates for the parameters are thus obtained by solving 

= (B {i) ,B obs ) , (4.6) 

3 

while the statistical error on the parameters is deduced from the second-order derivatives of 
X 2 , which define the Fisher matrix, given in our case by 

F« = (fl«,B tt ). (4.7) 

The Fisher matrix is a symmetric matrix, which can be determined by computing the 21 
different scalar products between the six elementary bispectra. 

For a real experiment, and if E-polarization is included as well, the above equations 
remain valid, except that the definition of the scalar product has to be replaced by a more 
complicated expression (see e.g. [45]): 



(5(0,^0)) s £ £ ^^(C^C^r^^BW^, (48) 



^■h loll 



where a, /3, 7, a ,/3 ,7 are polarization indices taking the two values T and E. The covariance 
matrix C; (a matrix in polarization space) is given by 

r - ( b i°T T + N F b i c F \ (a n 
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(CCC) 


(CCS) 


(C, SS) 


(s, CO 


(5,C5) 


(S,SS) 


3.9(2.5) x 10 -2 


4.5(3.6) x 1(T 2 
7.1(6.0) x 1(T 2 


2.3(2.1) x 1(T 4 
5.3(3.8) x 10~ 4 
28(6.4) x 10" 5 


2.4(1.6) x 10" 4 
3.8(2.1) x 10" 4 
16(3.7) x 1(T 5 
15(3.0) x 10~ 5 


6.9(4.3) x 10 -4 
11 (7.4) x 10~ 4 
33(9.5) x 10~ 5 
22(5.8) x 10~ 5 
5.1(1.6) x 10" 4 


5.3(3.1) x 10 -4 
8.8(5.5) x 10~ 4 
11(5.0) x 10" 5 
7.5(3.2) x 10~ 5 
2.4(1.0) x 10" 4 
21 (8.3) x 10" 5 



Table 1. Fisher matrix for the CDM isocurvature mode. Only the upper half coefficients are indicated, 
since the matrix is symmetric. The value between parentheses corresponds to the Fisher matrix 
components obtained without including the polarization. 



where bi is the beam function and Ni the noise power spectrum. We assumed the same 
beam function for temperature and polarization detectors, as well as no correlated noise, but 
the generalization is straightforward. In the calculation of the covariance matrix we only 
take the adiabatic power spectrum, since from observations we know that the isocurvature 
contribution to the power spectrum must be very small. 

For each type of isocurvature mode, we have computed the corresponding Fisher matrix 
by extending the numerical code described in [46] to include isocurvature modes and E- 
polarization, according to the expressions presented above. We have taken into account 
the noise characteristics of the Planck satellite [47], using only the 100, 143, and 217 GHz 
channels, combined in quadrature. Our computation goes up to Z max = 2500 and uses the 
WMAP-only 7- year best-fit cosmological parameters [13]. 

From the Fisher matrix, one can compute the statistical uncertainty on each of the 
parameters: 

Af = V / (^ rT k- (4-10) 

This takes into account the correlations between the various bispectra. By contrast, if one 
assumes that the data contain only a single elementary bispectrum, for example the purely 
adiabatic one, then the corresponding statistical error is 

A/ 1 = (single parameter) 

V Fa 

One can also determine the correlations between any two bispectra: 

c (^%- 

4.2 CDM isocurvature mode 

Our results for this mode have already been presented elsewhere [37], but we discuss here 
in more detail the peculiarities of the corresponding Fisher matrix, which is given in Ta- 
ble 1. One can immediately notice the intriguing fact that the coefficients of the upper left 
2x2 submatrix, corresponding to the purely adiabatic component and the correlated (£, (S) 
component, are typically two orders of magnitude larger than all the other coefficients. The 
correlation matrix, defined in (4.12) and given in Table 2 shows that the first two bispectra 
are strongly (anti-)correlated while their correlation with the four other bispectra is weak. 



(4.11) 



(4.12) 
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(CCC) (CCS) 


(C,ss) 


(£, CO 


(s,cs) 


(S, SS) 


1. -0.84(0.92) 
1. 


0.18(0.21) 
-0.15(0.18) 
1. 


0.003 (-0.12) 
-0.007 (+0.13) 
-0.069 (+0.24) 
1. 


-0.15(0.17) 
0.13(0.14) 
-0.80(0.97) 
-0.42 (0.43) 
1. 


0.13(0.18) 
-0.16(0.16) 
0.58(0.92) 
0.29(0.39) 
-0.82(0.98) 
1. 



Table 2. Correlation matrix for the CDM isocurvature mode. Only the upper half coefficients 
are indicated, since the matrix is symmetric. The value between parentheses corresponds to the 
correlations obtained without including the polarization (the absence of sign in the parentheses means 
that it is unchanged with respect to the value taking into account the polarization). 



From this Fisher matrix, one finds that the 68 % error on the parameters f l is given 

by 3 

Af = \J (F-^a = {9.6, 7.1, 160, 150, 180, 140} . (4.13) 

For ease of readability, we have written 160 instead of 1.6 • 10 2 , etc., but we are not claiming 
more than two digits of significance. We also remind the reader that in the purely adiabatic 
case, our f 1 , i.e. the (CCC) component, is 6/5 times the standard /nl- One sees that the 
first two uncertainties are typically one order of magnitude smaller than the last four. 

It is also interesting to estimate how much the inclusion of the polarization data in the 
analysis improves the precision of the non- linear parameters. The components of the Fisher 
matrix when the polarization is not taken into account can be read between the parentheses 
in Table 1. One notices that whereas the coefficients of the first two lines are reduced by a 
factor inferior to two, the other coefficients are significantly suppressed when one removes 
the polarization data. As a consequence, one finds that the uncertainties on the parameters 
without polarization, given by 

Af = \J (F~l)ii = {17, 11,980,390, 1060,700} (no polarization), (4.14) 

increase by less than a factor two for the first two parameters, whereas the increase is much 
bigger for the four other ones. 

The evolution of these uncertainties as a function of the cut-off Z max is shown in Fig. 15, 
both for the case where temperature and E-polarization data are used and for the case where 
only temperature data is included. One can see that the curves for temperature-only typi- 
cally look bumpier than the curves that include polarization as well. This can be explained 
as follows. First, unlike the power spectrum, the bispectrum is an alternating function, so 
that for certain regions in hhh space it is zero or close to zero, and the contribution to 
the determination of /nl> which is quadratic in the bispectrum, is then negligible in these 
regions. Second, as one can see for example in Fig. 12, the acoustic peaks of the polarization 
bispectrum are out of phase with the ones of the temperature one, so that including polar- 
ization neatly fills in the holes in Z1Z2Z3 space and leads to a smoother determination of /nLj 
as first pointed out by [48]. 

Our results can be understood by the following analysis in the squeezed limit (based on 
[46]), assuming that l\ = t <C I2 ~ ^3 = L. In this limit, the bispectra defined in (3.19) can 

3 The tiny differences in the 3rd, 5th, and 6th value compared to [37] are due to small improvements in the 
computer code. 
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Figure 15. Evolution of the /nl parameter uncertainties as one increases the cut-off l max , for the 
CDM isocurvature mode. The six thinner curves describe the situation if only temperature data is 
used, while for the six thicker curves both temperature and E-polarization data are included. 



be decomposed as 

POO 

bi'LL= / r 2 dr[a|(r)/3/(r)/3f(r) + ai(r)/3/(r)/3f(r) + ai(r)/3/(r)/3f(r)]. (4.15) 



The first term is subdominant, since, like the power spectrum, |/3/| decreases as 1/(1(1 + 1)) 
(or even faster as l/(l 2 (l + l) 2 ) for large I), as can be seen for example in Fig. 7. The last 
term, for instance, is explicitly given by 



2 X3 



7T 



kjdh k 2 dk 2 k 2 dk 3 g I L (k 1 )g J L (k 2 )gf(k 3 ) P((k 2 )P c (k 3 ) / r 2 dr j L (kir)j L (k 2 r)ji(k 3 r) , 



(4.16) 

where the last Bessel function oscillates slowly while the first two oscillate very rapidly. This 
leads to a cancellation of the radial integral unless k\ is very close to k 2 . We find that the 
above expression can thus be approximated by 

(f) 3 J kfdh J k\dk 2 J ^dfc 3 gi(fc 1 )gi(fc 2 )gf(fc3)P c (fc 2 )P c (A:3) 5 ^ 1 i ~ 2 fc2) ^ 00 yj e (k 3 r) 

J k 2 dkg[(k)gi(k)P c (k)^^ J k 2 dk 3 gf(k 3 )P Q (k 3 ) J™ yjifar) 

= Gu(L)U K (t)- ' (4.17) 

The S(k\ — k 2 ) is explained above, but together with the l/(kik 2 ) also motivated by the 
closure relation for spherical Bessel functions, r 2 jL(k\r)jL(k 2 r)dr = it8(k\ — k 2 )/(2kik 2 ). 
The 1/r follows from a dimensional analysis, and the 5 has been determined heuristically 
by comparing with the exact bispectrum: the ratio is approximately 5 and only weakly 
dependent on the small t. The full squeezed bispectrum (4.15) is thus approximated by 

t>iLL ~ Gu(L) H K (£) + g IK (L) Hj(£) . (4.18) 
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Figure 17. The absolute value of the squeezed bispectrum (multiplied by L(L + 1)) with l\ = 2 and 
l 2 = I3 = L plotted as a function of L. Both the exact bispectrum and (twice) the approximation 
given in (4.17) are plotted, for the purely adiabatic mode and the purely CDM isocurvature mode. 
The figure on the left shows the temperature-only TTT bispectra, the one on the right the pure 
polarization EEE bispectra. Note that the configuration is only squeezed in the case of large L; for 
small L the approximation should not be used. 



Assuming all primodial power spectra to be equal, the functions Qu are the angular (cross) 
power spectra plotted in Fig. 1 and 2. The function Hk(£) is simply an integral over /3f (r): 

f°° dr 

H K (l) = 5 -P?{r) (4.19) 
Jo r 

and is shown for small £ in Fig. 16. The good agreement of our approximation with the exact 
squeezed bispectrum is shown in Fig. 17. 

In the squeezed limit, one thus finds that only the first two elementary bispectra, (£, £C) 
and (CCS), depend on Q^(L). The four others depend on Q^s{L) and/or Qss{L). The 
large L limit of Q^s(L) and Q SS (L) are strongly suppressed with respect to Q^(L), which 
explains why the uncertainty on the first two non-Gaussianity parameters can be reduced 
by probing high multipoles (the bispectrum there is still sufficiently large compared to the 
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noise) while the uncertainty on the four other ones saturates as shown in Fig. 15. One 
can even understand why the curve for the (CCS) mode is below the one for (C ; CC) : their 
dominant terms both depend on the same Q^(L), but different %(£)-, and \Hs(&)\ > l^cWI- 
Finally, one can understand why including polarization helps much more for the uncertainty 
on e.g. the (S, SS) mode than for the ((, CO mode. As one can see from Fig. 15, it is in 
particular in the region 50 < / < 200 that the distance between the two (S, SS) curves 
increases compared to the distance between the two (£, CC) curves. A quick look at Fig. 1 
shows that in that region of multipole space the TT CDM isocurvature power spectrum (i.e. 
Qss{L)) becomes very small compared to the TT adiabatic spectrum (i.e. Q^(L)), but the 
EE CDM isocurvature spectrum still remains comparable to the EE adiabatic one. 
It is also instructive to compare (4.13) with the uncertainties 

Af = = {5.0, 3.7, 60, 82, 44, 69} (single parameter) (4.20) 

obtained by ignoring the correlations, or, equivalently, by assuming that only one parameter 
is nonzero. In particular, the contamination of the purely adiabatic signal by the other shapes 
induces an increase of the uncertainty, but only by a factor 2, which is rather moderate. 

Assuming that the adiabatic and isocurvature modes are uncorrelated implies that 
only the purely adiabatic and isocurvature bispectra are relevant. The corresponding two- 
parameter Fisher matrix, which is the submatrix of F^ with entries F\\, Fqq and Fiq, leads to 
uncertainties on and that are almost identical to the corresponding single-parameter 
errors. 

Finally, let us note that if the observed bispectrum is mainly purely isocurvature with 
amplitude f^ e \ a naive analysis using only the purely adiabatic estimator would lead to an 
apparent adiabatic coefficient 

/ ~(l) = ^6 / ~(6)^ 10 -2 / ~(6) ; (421) 

Fu 

thus hiding the isocurvature signal with larger amplitude. 
4.3 Baryon isocurvature mode 

The Fisher matrix for the baryon isocurvature mode can be easily deduced from the CDM 
Fisher matrix. Indeed, as discussed earlier, the CDM and baryon isocurvature transfer func- 
tions are identical up to a rescaling by u;& c introduced in (2.11). Consequently, the a and /3 
functions are simply rescaled: 

af b (r) = u bc af»(r), tf h (r) = uj hc p^ (r) . (4.22) 

The rescaling of the various bispectra will thus depend on the number of S indices, i.e. 

bi;^ 3 (s b ) = (u bc )ni;ll(s c ), (4.23) 

where the power p is the number of S among the indices {UK}. In summary, all coefficients 
of the baryon isocurvature Fisher matrix can be deduced from Table 1 by using the rescaling 

F§> = Mi MjF?f (no summation), M t = {l,u bc , uo 2 bc , u bc , uo 2 bc , wfj, (4.24) 

where, in our computation, u) bc = 0.2036. 
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(CCC) 


(CCS) 


(C, SS) 


(s, CC) 


(S,(S) 


(S,SS) 




3.9(2.5) x 10 -2 


3.6(2.6) x 1(T 2 


5.6(4.1) x 10 -3 


7.9(5.1) x 10 _a 


8.8(6.1) x 10 -3 


2.2(1.6) x 10" 


-3 




3.8(2.9) x 1(T 2 


6.3(4.8) x 10" 3 


7.4(5.2) x 1CT 3 


9.2(6.7) x 10" 3 


2.5(1.8) x 10" 


-3 






11 (8.1) x 1CT 4 


12(8.5) x 10"" 4 


1.6(1.1) x 1(T 3 


4.4(3.1) x 10" 


-4 








1.8(1.1) x 10" 3 


2.0(1.3) x 10"~ 3 


5.0(3.2) x 10" 


-4 










2.5(1.6) x 10" 3 


6.8(4.4) x 10" 


-4 












2.1 (1.2) x 10" 


-4 



Table 3. Fisher matrix for the neutrino density isocurvature mode. Only the upper half coefficients 
are indicated, since the matrix is symmetric. The value between parentheses corresponds to the Fisher 
matrix components obtained without including the polarization. 



The parameter uncertainties can also be deduced from the CDM results via the rescal- 
ings Aft: Af(S b ) = A/ i (5 c )/A/' i . One thus obtains: 

Af = {9.6, 35, 4000, 720, 4300, 16600} . (4.25) 

Except for the purely adiabatic coefficient, we thus find that the uncertainties on all the other 
coefficients are significantly larger than the uncertainties obtained in (4.13) in the CDM 
case, simply because the elementary bispectra have a smaller amplitude than their CDM 
counterparts. By contrast, the correlation matrix, which is independent of the normalization 
of the bispectra, is exactly the same as in the CDM case. 

4.4 Neutrino density isocurvature mode 

For a neutrino density isocurvature mode, we have obtained the Fisher matrix in Table 3. 
Unlike the case of CDM isocurvature, here the difference between the different entries in 
the Fisher matrix is smaller, although the coefficients in the upper left 2x2 submatrix are 
still about one order of magnitude larger than the others. Also in contrast to the CDM 
isocurvature case, we see that all coefficients increase about equally when polarization is 
included. 

The corresponding uncertainties on the six non-Gaussianity parameters are (taking into 
account the correlations) 4 

Af = {28, 36, 190, 150, 240, 320}. (4.26) 

When using temperature only, the uncertainties increase to 

Af = {58, 75, 540, 340, 720, 950} (no polarization). (4.27) 

The evolution of the uncertainties as a function of Z max is shown in Fig. 18. As in the 
CDM case, the (£, CC) an d (C> £S) non-Gaussianity parameters can be determined more 
accurately than the other four, although the difference is not as big as for CDM. Unlike 
for CDM, all parameters gain about the same from the inclusion of polarization, and all 
uncertainties continue to decrease when higher multipoles are probed, since the neutrino 
density isocurvature power spectrum does not decrease as steeply as the CDM isocurvature 
one. 

4 While we were finalizing our manuscript, we became aware of the work [49], where the authors also 
investigate neutrino density isocurvature non-Gaussianity. Their numbers for the uncertainties are very similar 
to ours (note that they use the six non-Gaussianity parameters that we introduced in [37] but in a different 
ordering), although they use a different selection of Planck channels. 
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(CCC) (CCS) 


(C,ss) 


(S, CO 


(s, cs) 


(S, SS) 


1. -0.83(0.82) 
1. 


0.42 (0.31) 
-0.58(0.45) 
1. 


-0.73(0.67) 
0.58 (0.55) 
0.03 (0.35) 
1. 


0.26(0.17) 
-0.31 (0.22) 
-0.45 (0.69) 
-0.76 (0.81) 
1. 


-0.16(0.09) 
0.25 (0.20) 
0.19(0.41) 
0.58 (0.64) 

-0.84 (0.86) 
1. 



Table 4. Correlation matrix for the neutrino density isocurvature mode. Only the upper half coef- 
ficients are indicated, since the matrix is symmetric. The value between parentheses corresponds to 
the correlations obtained without including the polarization (the absence of sign in the parentheses 
means that it is unchanged with respect to the value taking into account the polarization). 
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Figure 18. Evolution of the /nl parameter uncertainties as one increases the cut-off Z max , for the 
neutrino density isocurvature mode. The six thinner curves describe the situation if only temperature 
data is used, while for the six thicker curves both temperature and E-polarization data are included. 

The correlation matrix is given in Table 4. If one assumes the parameters to be inde- 
pendent, one finds 

Af = {5.0, 5.1, 30, 24, 20, 69} (single parameter). (4.28) 

One sees that the increase of the uncertainties due to the correlations is much more important 
here than for CDM, due to the larger correlations between the various modes. 

4.5 Neutrino velocity isocurvature mode 

For a neutrino velocity isocurvature mode, we have obtained the Fisher matrix in Table 5. 
One notices that, including polarization, all entries are roughly of the same order of magni- 
tude, but without polarization, they vary a lot. The corresponding uncertainties on the six 
non-Gaussianity parameters are (taking into account the correlations) 

Af = {25,22,85,81,77,71}. (4.29) 
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(CCC) 


(CCS) 


(C,ss) 


(SCO 


(SCS) 


(S SS) 


3.9(2.5) x 10 -2 


3.5(1.1) x 1(T 2 


74(5.1) x 10 -4 


15(9.5) x 10" 3 


16(5.6) x 10" 3 


40(2.6) x 10" 4 




73(7.6) x 1(T 3 


23(1.3) x 10~ 3 


13(4.0) x 10~ 3 


30(3.4) x 10~ 3 


12(0.65) x 10~ 3 






80(4.3) x 1(T 4 


29(1.6) x 10" 4 


96(4.8) x 10" 4 


45(2.1) x 10~ 4 








5.9(3.6) x 10~ 3 


6.3(2.1) x 10" 3 


16(0.88) x 10~ 4 










13(1.6) x 10" 3 


54(2.5) x 10" 4 












28(1.1) x 10" 4 



Table 5. Fisher matrix for the neutrino velocity isocurvature mode. Only the upper half coefficients 
are indicated, since the matrix is symmetric. The value between parentheses corresponds to the Fisher 
matrix components obtained without including the polarization. 



When using temperature only, the uncertainties increase to 

Af = {51, 120, 460, 180, 190, 570} (no polarization). (4.30) 

The evolution of the uncertainties as a function of Z max is shown in Fig. 19. 

One sees that in this case the difference between the first two and the other four un- 
certainties (when including polarization) is much smaller than for CDM or neutrino density 
isocurvature (a factor of about 3 compared to a factor of about 15 in the CDM case). In 
particular, the latter four can be determined more accurately than in the case of CDM or 
neutrino density isocurvature. However, the improvement due to polarization is much more 
important than in the neutrino density case, in particular for the (£, (S), the ((,SS) and 
the (S, SS) parameters. One can understand why the (£, CS) mode, for example, gains much 
more from polarization than the (£, ££) mode with a similar argument as the one presented 
for CDM isocurvature. The dominant contributions to these modes both depend on Q^(L) 
(defined in (4.17)), but for (C, CO this is multiplied by U c (i) and for (£ (S) by Hs(£). And as 
one can see from Fig. 16, the ratio Hs/Hc; for neutrino velocity isocurvature increases enor- 
mously when one passes from temperature to polarization (remember that it is the lowest 
values of t that contribute most to the squeezed configuration) . 

The correlation matrix is given in Table 6. If one assumes the parameters to be inde- 
pendent, one finds 

Af = {5.0, 3.7, 11, 13, 8.7, 19} (single parameter). (4.31) 

Hence one sees that the correlations in the case of neutrino velocity isocurvature are more 
important than for CDM, but less than for neutrino density. 

5 Constraints on early universe models 

In the previous section, we have studied how to obtain constraints on the six non-linearity 
coefficient f NL without assuming any particular relation between them. In the context of an 
early Universe model, or in a class of models, one can go further and use the results of the 
previous section to obtain some constraints on the parameters of the model itself. 

5.1 General analysis 

As we have seen earlier, isocurvature perturbations require the existence of at least two 
degrees of freedom in the early Universe. So, for simplicity, let us focus on models with 
two scalar fields, (j) and a, such that isocurvature perturbations are generated only by the 
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(CCC) (CCS) 


(C,ss) 


(s, CO 


(S,(S) 


(S, SS) 


1. -0.50(0.42) 
1. 


0.08 (0.11) 
-0.60(0.80) 
1. 


-0.82(0.66) 
0.25 (-0.18) 
0.37(0.44) 
1. 


0.34 (0.34) 
-0.26(0.35) 
-0.51 (+0.07) 
-0.74(0.61) 
1. 


-0.13 (+0.05) 
0.36 (0.39) 
-0.29(0.75) 
0.25 (-0.10) 
-0.46(0.36) 
1. 



Table 6. Correlation matrix for the neutrino velocity isocurvature mode. Only the upper half 
coefficients are indicated, since the matrix is symmetric. The value between parentheses corresponds 
to the correlations obtained without including the polarization (the absence of sign in the parentheses 
means that it is unchanged with respect to the value taking into account the polarization). 



Errors on fNL for neutrino velocity isocurvature 




Figure 19. Evolution of the /nl parameter uncertainties as one increases the cut-off Z max , for the 
neutrino velocity isocurvature mode. The six thinner curves describe the situation if only temperature 
data is used, while for the six thicker curves both temperature and E-polarization data are included. 

fluctuations of a and all non-linearities are also dominated by their contribution from a. This 
means that we have 

C = Nfa + N<5a + l -N<Ja 2 , S = N s a 5a + ~iV£> 2 . (5.1) 

Using the general expressions (3.15-3.17), one easily finds that the coefficients are in- 
terdependent and can be written in the form 

S£ C = M/3 2 , (5-2) 

/£f = £<s Ml« 1/2 S 3 / 2 , (5.3) 

/jVi (5.4) 
where we have introduced the contribution of a in the adiabatic power spectrum, 

(ivS) 2 



(iVi) 2 + (Nfr 



(5.5) 
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and the isocurvature to adiabatic ratio, 



a = } N ° )2 t . (5.6) 
(JVj) 2 + (JV<)= 

Note that the extraction of an isocurvature component in the power spectrum would fix a 
while, in principle, H could be determined from observations by measuring both the bispec- 
trum and the trispectrum coefficients, since they satisfy consistency relations [36] similar to 
the purely adiabatic relation tnl = f% L /E. 
The two coefficients 

M/ = ^V" (I = C,S) (5.7) 

fully characterize the non-Gaussianity of the adiabatic and isocurvature perturbations, re- 
spectively, while e fs denotes the relative sign of C and S: e cs = 1 if they have the same 

sign, e cs = —1 otherwise. Interestingly, /jy£ and /j^f 5 share the same sign as fij, whereas 

flfi can have a different sign. The hierarchy between the coefficients fl}^ ' , I being fixed, 

~I CC ~I ss 

depends on the relative amplitude of a and S: /jvl; dominates if H S> a, whereas f^ L 

dominates if E <C a. 

For given values of a and E, the uncertainties on the two parameters ^ and fig are 
determined from the "projected" Fisher matrix 



with 

7r Wf = { ~ 2 , e^a 1 ^ 2 , aE, 0, 0, 0}, n {i)s = {0, 0, 0, H 2 , e^a 1 ' 2 ??' 2 , aE} . (5.9) 

From this 2x2 Fisher matrix, one can easily deduce the expected uncertainties on the two 
parameters ^ and /is, by using the analog of (4.10). 

5.2 Illustrative example 

In [35, 37] we have studied a class of models which produces perturbations of the above 
type. In these models, a is a curvaton which decays into radiation and CDM. Since part of 
the CDM can have been produced before the decay of the curvaton, one can introduce as a 
parameter the fraction of CDM created by the decay as 

/c= 7 ;^ , (5.io) 

i l c + 7 c i No- 
where the O's represent the relative abundances just before the decay and j c is the fraction 
of the curvaton energy density transferred into CDM. The second relevant parameter, 

-- V -■*>"• (5.11) 



(l- 7c S!„)(4-fl„) 

quantifies the transfer between the pre-decay and post-decay perturbations [35] (one finds 

.after = (1 _ r)C bc 



Rafter = (1 _ r ) ^before &t the linear leyel ) 
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As shown in [35, 37], one can derive the "primordial" perturbations £ and S as expan- 
sions, up to second order, in terms of 5a, which yield the coefficients and N* a . Using 
these results and assuming r<l, one obtains 

K = Y r > MS = ^2(/ c (l-2/ c )-r), (5.12) 
while e (S = sgn(/ c — r) and 




(5.13) 



In the regime / c < r < 1, one finds = 3/(2r), ^5 = — 9/(2r) and a = 9H, with e cs = — 1 
and the amplitudes of non-Gaussianities depend only on the parameter r. By contrast, in 
the opposite regime r <C / c *C 1, ^ = 3/(2r), /j*s = 9/ c /(2r 2 ) and a = 9E(f c /r) 2 . The 
coefficients f N ^ J , which depend on [Ag, are thus enhanced with respect to the coefficients 
f^jf in the latter case. 

6 Conclusion 

We have systematically investigated the angular bispectra generated by initial conditions 
that combine the usual adiabatic mode with an isocurvature mode, assuming local non- 
Gaussianity. We have studied successively the four types of isocurvature modes, namely cold 
dark matter (CDM), baryon, neutrino density and neutrino velocity isocurvature modes. In 
each case, the total bispectrum can be decomposed into six elementary bispectra and we have 
estimated the expected uncertainties on the corresponding coefficients, which are extensions 
of the usual purely adiabatic /nl parameter, in the context of the forthcoming Planck data 
of the cosmic microwave background radiation (CMB). As we showed, the results for baryon 
isocurvature can be obtained from a simple rescaling of the CDM isocurvature results, but 
the others are distinct. 

In the squeezed limit, where one multipole I is much smaller than the other two (which 
then have to be almost equal due to the triangle inequality), we have shown that the six 
elementary bispectra factorize as a function of the small I times the power spectrum as a 
function of the large I. Since the squeezed limit components dominate the bispectrum for 
local non-Gaussianity, we have been able, using this factorization, to give simple explanations 
for the various interesting results that we have observed. 

By enlarging the space of initial conditions, one obviously expects a larger uncertainty 
on the purely adiabatic /nl coefficient. Interestingly, this uncertainty is increased only by 
a factor 2 for CDM and baryon isocurvature modes, whereas it increases by a factor 6 in 
the neutrino density isocurvature case and by a factor 5 in the neutrino velocity isocurvature 
case. This can be explained by the fact that the CDM isocurvature power spectrum decreases 
much faster with I than the adiabatic and neutrino isocurvature ones. As we have shown, this 
means that the uncertainties on jp'" and f^^ s in the case of CDM isocurvature continue to 
improve as one increases the number of available multipoles, while the other four saturate 
at a much lower I. As a consequence the first two can be determined much more accurately 
than the other four, and are only weakly correlated with them. This small correlation also 
means that it is important to look at the data with the full /nl estimator, and not just the 
adiabatic one, as a large CDM isocurvature non-Gaussianity could be hiding behind a small 
adiabatic signal. 
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We have shown that the E-polarization often plays a crucial role in reducing the uncer- 
tainties. In the CDM isocurvature case, polarization improves slightly the precision on the 
coefficients and f^'^ s , but the precision of the other four coefficients improves by a factor 
of order five. Polarization is also very important for some of the parameters in the neutrino 
velocity case, the uncertainty on the purely isocurvature f s,ss , for example, improving by a 
factor 8 when polarization is included. Again we were able to explain these results from the 
behaviour of the power spectra, using the factorization of the squeezed bispectrum. 

The decomposition of the bispectrum into six elementary bispectra and the CMB con- 
straints on the six parameters do not depend on any assumptions about the specific 
model of the early Universe. We only assumed that (possibly correlated) primordial adia- 
batic and isocurvature modes are produced, with a primordial bispectrum of local type and 
with power spectra that all have the same shape. Note that neither of these assumptions 
appears to be essential; they were only made for simplicity. In future work we will look into 
the possibility of relaxing these assumptions. 

If, however, one does consider an explicit early Universe model, there are often relations 
between the six parameters, and an observational detection or constraint could then be used 
to check for such a relation and put constraints on the parameters of the model. We have 
discussed a general class of models with two scalar fields, where only one of the fields generates 
both the isocurvature perturbations and all non-linearities. We have also considered a specific 
implementation of this general model, where the two fields are an inflaton and a curvaton. In 
this model, a CDM isocurvature mode is produced and the six parameters only depend 
on two model parameters. In some ranges of the model parameters, the isocurvature mode 
is subdominant in the power spectrum but provides observable non-Gaussianity that can 
dominate the usual adiabatic non-Gaussianity. Looking for these new angular shapes in the 
CMB data would thus provide interesting information on the very early Universe. 
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